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The onset of convection in the form of magneto-inertial waves in a rotating fluid sphere 
permeated by a constant axial electric current is studied through a perturbation anal- 
ysis. Explicit expressions for the dependence of the Rayleigh number on the azimuthal 
wavenumber are derived in the limit of high thermal diffusivity. Results for the cases 
of thermally infinitely conducting and of nearly thermally insulating boundaries are ob- 
tained. 



1. Introduction 

Buoyancy driven motions of an electrically conducting fluid in a rotating system and 
in the presence of a magnetic field represent a fundamental aspect of the dynamics in 
stars and in planetary interiors. It is thus not surprising that this topic has received 
much attention ever since Chandrasekhar (1961) and others have presented the first 
mathematical solutions for problems of this kind. Because of the numerous parameters 
entering these problems it is not always easy to gain an overview of the physically relevant 
convection modes. Analytical expressions are thus especially desirable, even though they 
can usually be obtained only for the onset of convection in planar layers. For some 
examples of this kind we refer to Chandrasekhar's monograph (1961). From the geo- and 
astrophysical point of view, magnetoconvection in spherical configurations is of special 
interest. Here numerical schemes are usually required for the solution of the mathematical 
problems. For an overview of the subject we refer to the article by Fearn (1994). 

In the case of low Prandtl numbers convection typically assumes the form of inertial 
waves. Buoyancy and viscous dissipation can be treated as perturbations and analytical 
expressions for the critical values of the Rayleigh numbers for onset of convection can 
be obtained. This approach has been employed in the case of a rotating sphere without 
magnetic field first by Zhang (1994). For later analyses with different boundary conditions 
see Zhang (1995) and Busse & Simitev (2004). 

Magneto-inertial modes can be excited in the presence of both, the Coriolis force and 
the Lorentz force of an imposed magnetic field. Of particular interest from an astro- 
physical point of view are magneto-inertial waves in rotating fiuid spheres with a strong 
azimuthal magnetic field. Malkus (1967) has demonstrated that for a magnetic field cor- 
responding to an axial electric current with constant density the dispersion relation for 
magneto-inertial modes can be expressed in terms of the dispersion relation of inertial 
modes in the absence of a magnetic field. He thus obtained simple analytical expressions 
for magneto-inertial waves in a configuration of special interest for geophysical and as- 
trophysical applications. In this paper we use the Malkus approach as a starting point 
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Figure 1 . Geometrical configuration of tlie problem. A part of the outer spherical surface is 
removed to expose the interior of sphere to which the conducting fluid is confined. 

for an analysis of the onset of convection in a rotating sphere permeated by a azimuthal 
magnetic field. In analogy to the treatment of inertial convection in the paper of Busse 
and Simitev (2004) we investigate the effect of buoyancy and ohmic dissipation as pertur- 
bations on magneto-inertial waves. The main result is the determination of expressions 
for the critical Rayleigh number for the onset of convection as a function of the rotation 
rate, the strength of the magnetic field and the other parameters of the problem. 

In the following we start with the mathematical formulation of the problem in section 
2. The special limit of a high ratio of thermal to magnetic diffusivity will be treated 
in section 3. The general case requires the symbolic evaluation of lengthy analytical 
expressions and will be presented in section 4. A discussion of the results and an outlook 
on related problems will be given in the final section 5 of the paper. 

2. Mathematical formulation of the problem 

We consider a homogeneously heated, self-gravitating fluid sphere as shown in figure 
[TJ The sphere is rotating with a constant angular velocity flk where k is the axial 
unit vector. The fluid is electrically conducting and is characterized by its magnetic 
diffusivity 77, its kinematic viscosity ly and its thermal diffusivity k. Following Malkus 
(1967) we assume that the fluid sphere is permeated by a magnetic field of the form 
B = Bok X r where r is the position vector with respect to the center of the sphere and 
r is its length measured in fractions of the radius tq of the sphere. Since the Lorentz 
force like the centrifugal force can be balanced by the pressure gradient a static state 
exists with the temperature distribution Ts = Tq — Pr^r^ /2 and the gravity field given 
by g = —^r^r. We use the length tq, the time r-g/yy, the temperature rf/^arQ and 
the magnetic flux density ^JJ^QTj/rQ as scales for the dimensionless description of the 
problem. The density is assumed to be constant except in the gravity term where its 
temperature dependence given by a = {dg/dT)/g = const, is taken into account. Since 
we are interested in the onset of convection in the form of small disturbances we can 
neglect all quadratic terms in the equations. The linearized versions of the equations of 
motion, the equation of magnetic induction, and the heat equation for the deviation & 
from the static temperature distribution are thus given by 



dtu + Tkxu + V{n - b ■ j x r) + {j x r) ■ Vb - j x b ^ Qr + P^V^it, (2.1&) 



V-it = 0, V-b = 0, 



(2.1a) 



dtb — {j X r) ■ Vii + j x u — V^b, 



(2.1c) 
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Figure 2. Competition of modes with increasing tP in the non-magnetic case discussed in Busse 
& Simitev (2004). The Rayleigh number R as a function of tP for m = (thick dash-dotted 
hnes) and m = 1 (thin lines). Results based on the explicit expressions (4.6) and (3.4) from 
Busse & Simitev (2004) are shown in solid lines and broken lines respectively in the case m = 1. 
(a) Case A, fixed temperature boundary conditions, (b) Case B, insulating thermal boundary 
conditions. 



Rr u + V^e - PdtQ/Prn = 0, 



(2. Id) 



where the dimensionlcss magnetic field has been assumed in the form j x r + b. Here 
2j = 2jk is the vector of the density of the imposed electric current. The Rayleigh number 
R, the Coriolis parameter r and the Prandtl number P and the magnetic Prandtl number 
Pm are defined by 



R 



2nr 



P ~ ~l Pm ~ ~, S = — = — . 

K Tj K _r,fj 



(2.2) 



We have attached a " to the Rayleigh number to remind the reader that the magnetic 
diffusivity replaces the kinematic viscosity in the definition of R. We have also introduced 
the Strebor number S which is the inverse of what often is called the Roberts number. 
Without losing generality we may assume an exponential dependence on time t and on 
the azimuthal angle (/>. Since both, the velocity field u and the magnetic field b, are 
solenoidal we may use the the representation of these vector fields in terms of poloidal 
and toroidal components. 



il = tXCXp (i{ljJTt + 7710)) 

b ~ bexp (i{ujTt + rn^)) 



Equation (|2.1cp can now be written in the form iurb 
eliminate b from equation (12. 16^ . 



= (V X (Vw X r) + Vw X r) exp [i(ujTt + mcj))), (2.3a) 
: X (V/7 X T*) + V(7 X T*) exp {i{ujTt + m(j>)) , (2.36) 
e = e exp {i{LUTt + mc/))) . (2.3c) 

V^6 which allows us to 



imju 



7 ] k X u — Vtt 

-t-— ^fc X V^it 



1 



-Or 



? 9 



— fc X V^bfc 



P„ 



(2.4) 



ILOT T 

derivative is replaced by im, of 



where 7 is defined by 7 = j/t. In the V-operator the 
course. In equation (|2.4p the magnetic field b appears only in the form of the boundary 
layer correction 6b which is required since the basic dissipationless solution does not 
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satisfy all boundary conditions. For the same reason the Ekman layer correction Ub had 
to be introduced. 

In the following we shall assume the limit of large r such that in first order of approxi- 
mation the right hand side of equation (|2.4p can be neglected. The left hand side together 
with the condition V-it — is of the same form as the equation for inertial modes. Among 
the latter those corresponding to the sectorial spherical harmonics usually lead to the 
lowest critical Rayleigh numbers for the onset of convection. We shall assume that this 
property continues to hold as long as 7^ is sufficiently small. The sectorial inertial modes 
are given by 



vo = P:^ {cos 0)f{r), iz;o = Pr+i(cos0)V(r), (2.5) 



with 



/(r) =r™-r"+2, ^(r)=r"+i^ 2.m(m + 2) ^ ^ 

' v-V ; (2m + l)(Ao(m2 + 3m + 2)-m)' ^ ' 



2^ + 3 J' ^'-''^ 

where Aqt is the frequency of the inertial modes. The sectorial magneto-inertial modes 
are described by the same velocity field (|2.5p and their magnetic field is given by bo = 
mjUQ/ujQ. As before the subscript refers to the dissipationless solution of equations 
(|2.1b .b.c). The frequency luq of the magneto-inertial waves is determined by 

9 99 

Ao = (2.7) 



which yields 



UJQ 



Aq , / X? 



2 ±y^+m72(m-Ao). (2.8) 



With account of (2.6b), this dispersion relation allows for a total of four different fre- 
quencies Wo- For small values of 7^ these are given by 



1 / , , + 4to + 3 \ 2 2 I 1 I / + 4to + 3 1 2 



-1 



The upper sign in expression (|2.9ap refers to retrogradely propagating modified inertial 
waves, while the lower sign corresponds to the progradely traveling variety. The effect of 
the magnetic field tends to increase the absolute value of the frequency in both cases. 
Expression (|2.96|) describes the dispersion of the slow magnetic waves. The upper sign 
refers to the progradely traveling modified Alfven waves and the lower sign corresponds 
to retrogradely propagating modified Alfven waves. 

The magneto-inertial waves described by expressions (|2.5p satisfy the condition that 
the normal component of the velocity field vanishes at the boundary. This property 
implies that the normal component of the magnetic field vanishes there as well. Additional 
boundary conditions must be specified when the full dissipative problem described by 
(|2.4p is considered. We shall assume a stress- free boundary with either a fixed temperature 
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(case A) or a thermally insulating boundary (case B), 

, , , ,16 = case A 

r • M = r • V(r X «)/?'= and < at r = 1. (2.10) 

I drO = case B 

Additionally we shall assume an electrically insulating exterior of the sphere which re- 
quires 

g = at r = l (2.11) 

and the matching of the poloidal magnetic field to a potential field outside the sphere. 

Following the procedure of earlier papers (Zhang, 1994; Bussc and Simitev, 2004) we 
use a perturbation approach for solving equations (j2.ip in the limit of large r, 



u~Uq + Ui + ..., uj ~ ujo + oji + ... (2-12) 

After the ansatz ()2.12|) has been inserted into equation ()2.4|) such that terms with ui 
appear on the left hand side, while those with uq and toi appear on the right hand side, 
we obtain the solvability condition for the equation for Ui by multiplying it with Uq and 
averaging it over the fluid sphere, 

= (Or . n*) + ((nS • V^«o)^ + («S • V^6o.)) - 4^^) \ 

where the brackets (...) indicate the average over the fluid sphere and the * indicates the 
complex conjugate. We have neglected all terms connected with viscous dissipation, i.e. 
we have assumed the of vanishing P^, since we wish to focus on the effect of ohmic dissi- 
pation. The effects of viscous dissipation have been dealt with in the earlier paper (Busse 



and Simitev, 2004). Since (mJ • V ito) vanishes, as demonstrated in Zhang ct al. (2001) 



we must consider only the influence of the boundary layer magnetic field bob- It is deter- 
mined by the equation 

iujoTbob = y^bob. (2.14) 
Since the solutions of this equation are characterized by gradients of the order -i/r, the 
boundary layer correction needed for the poloidal component is of the order -y/r smaller 
than the correction needed for the toroidal component. For large r we need to take into 
account only the contribution got given by 



.go(r = l)exp (-(1 + zs)(l - r)/^^) 

wo{r = 1) exp (-(1 + is)(l - ?Vko|T/2) , (2.15) 



9ob = - 

where s denotes the sign of loq. The solvability condition thus becomes reduced to 



T 2 (oJo - 77172)y/2|wo|t J-1 



X (m + 1)(to 



2m(m + 2) 



(2'^ + 1) (^^^("^ + l)(m + 2) - m) 
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Figure 3. Case A. (a) The critical Rayleigh number Rc as a function of the wave number m for 
7 = 0.1 and r = 10^ . . . 10® increasing from bottom with log-scale decades given by thick lines, 
(b) The critical Rayleigh number Rc as a function of 7 for S — 1 and m = 1 . . .8 increasing from 
bottom, (c) Competition of modes with increasing S for 7 — 0.1 and m = 1. Explicit expressions 
1)3. 20p in the limit tS <^ 1 are shown by broken lines, (d) The critical Rayleigh number Rc as 
a function of S for 7 = 0.1 and m = 1 ... 8 increasing from bottom. The axisymmetric mode 
m = is given for comparison in panels (c,d) by a dot-dashed line. In panels (b) to (d) r = 10*. 



3. Explicit expressions in the limit tS <^ 1 

The equation (|2.1(i[l for O can most easily be solved in the limit of vanishing loqtS. In 
this limit we obtain for O, 

e = {cos 9) exp{irrup + iujTt)q{r), (3.17) 

with 

-/ m(m + iy+-^ m(m + l)r'"+2 

q(r) = R I — — CT 

+ 5){m + 'i) - {m + l)m {m + S){m + 2) - {m + l)m 

(3.18) 

where the coefficient c is given by 

1 1 

case A, 



(m + 5)(m + 4) - (m + l)m (m -|- 3)(m + 2) - {m + l)m 

, , , , . , ^3.1Q) 

(m + 4) /m [m -\-2) m 

- — - CclSG o 

(m + 5)(m -f 4) - (m -I- l)7Ti {m + i){m + 2) - [m + l)m 
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Since Q and the left hand side of equation (2.19) is imaginary, the real parts of the two 
terms on the right hand side must balance. Wc thus obtain for R the result 



i? = s 



r 7^(m-a;o) 
2\u!a\ {uo - mj'^) 



m{m + 2) 



-(m + l)(m + 2) — m 



X (2m + 9)(2m + 7)(2m + 5)^{2m + 3) 



(3.20) 
m + 2 1 



where the coefficient b assumes the values 



{m{10m + 27) case A, 

14m2 + 59™ + 63 case B. 



771+ 1 6' 



(3.21) 



Obviously the lowest value of R is usually reached for m ~ 1, but the fact that there 
are four different possible values of the frequency complicate the determination of the 
critical value Rc. Expression (|3.20p is also of interest, however, in the case of spherical 
fluid shells when the (777 — l)-modc is affected most strongly by the presence of the inner 
boundary. Convection modes corresponding to higher values of 777 may then become 
preferred at onset since their r-dependcncc decays more rapidly with distance from the 
outer boundary according to relationships (|2.6p . 
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4. Solution of the heat equation in the general case 

For the solution of equation ()2.1d|) in the general ease it is eonvenient to use the Green's 
function method. The Green's function G(r, a) is obtained as solution of the equation 

[dyOr + ( - iuJorS - m{m + 1))] G(r, a) = 6{r - a), (4.22) 

which can be solved in terms of the spherical Bessel functions jmip-r) and y„i{fir), 

Gi{r, a) = Aijmil^r) for ^ r < a, 

G2{r, a) = Aj,n{nr) + By^ifir) for a < r ^ 1, 



G(r,a) 



(4.23) 



where 



A = -/ij,„(Afa) 



B = ^-ijmip-a). 



A solution of the equation (j2.1(ij) can be obtained in the form 

r-l 



(4.24a, b) 
(4.24c, d) 



g(r) = — m(m + l) / G2{r,a)[a 



Evaluations of these integrals for to = 1 yield the expressions 



)a2da+ / Gi(r,a) (a"-a"+2^a2da 




r(/x2 + 10) - fi^r^ - 



10(/ir cos(/ir) — sin(^r)) \ 
r2(/iCOS/i — sin/i) / 

(/i^ — 10) (/^r cos(/xr) — sin(/xr)) \ 



(4.25) 



case A, 



case B. 



r(/i2 + 10) - /iV .^ • A 

[2fi cos /i — (2 — /i^) sm /ij 

(4.26) 

Lengthier expressions are obtained for to > 1. Expressions (|4.26p can now be used to 
calculate R and uji on the basis of equation (|2.16p . In the case to = 1 we obtain 

/2^72 (c^o - 1) 



189 
"20" 



^\uj^\ {ujo - 7^) (6 Ao - 1)^ 



/i"^ - 525^4"** - 175 Re 



sm jj, 



/i^(/i cos fi — sin ji) 



X < 



231/i" 



7Re. 



(/i^ — S^-^ + 9^) cos ^ — 9 sin ^ 
fj.^[{fj.'^ — 2) sin^ + 2^cos^) 



(4.27) 



case A, 



case B, 



where Re{} indicates the real part of the term enclosed by {}. Expressions (|4.27p have 
been plotted as functions of S in figures[3fc) and[4{c) for the cases A and B, respectively. 
Four curves appear since there are four possible values of ujq for each to. For values S 
of the order 10~^ or less, expressions p.20p are well approached. The retrograde mode 
corresponding to the positive sign in (j2.6b|) always yields the lower value of R but it 
looses its preference to the progradely traveling modified Alfven mode corresponding to 
the upper sign in (j2.96p as S becomes of the order 10~^ or larger. This transition can 
be understood on the basis of the increasing difference in phase between 6 and Ur with 
increasing S. While the mode with the largest absolute value of uj is preferred as long 
as Q and Ur are in phase, the mode with the minimum absolute value of uj becomes 
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Figure 5. The border where the transition from modes characterised by ujo\ to modes charac- 
terised by cJo3 occurs in various sections of the parameter space. The value of the parameters 
are m = 1, S = 1, 7 = 0.1, and r = 5000 where they are not varied on the axes. Case A is 
denoted by a solid lines and Case B by broken lines. 



preferred as the phase difference increases since the latter is detrimental to the work 
done by the buoyancy force. The frequency perturbation uji usually makes only a small 
contribution to lo which tends to decrease the absolute value of lo. This transition shifts 
towards smaller values of S and 7 as r is increased as illustrated in figurcEl The magneto- 
inertial convective modes corresponding to higher values of m = 1 ... 8 exhibit similar 
behaviour as figures [31[d) andUJd) demonstrate for the cases A and B, respectively. The 
value m = 1 is always the preferred value of the wave number, except possibly in a very 
narrow range near 7 = 0.03 as indicated by figure [nta,b) in the case A and possibly 
near 7 = 0.02 in the case B and figure |31Ia,b). The axisymmetric mode m = 0, given 
for comparison in panels (c) and (d) in figures [3] and 31 is never preferred in contrast to 
the purely non-magnetic case where it becomes the critical one near the transition from 
retrograde to prograde inertial convection modes as seen in figure [2l 

For very large values of t and S the Rayleigh number R increases in proportion to 
y/rirSY for fixed m. In spite of this strong increase remains of the order t^^'^S on 
the right hand side of equation (j2.1bp . The perturbation approach thus continues to be 
valid for r — > 00 as long as 5 <C 1 can be assumed. For any fixed low value of 5", 
however, the onset of convection in the form of prograde inertial modes will be replaced 
with increasing t at some point by the onset in the form of columnar magneto-convection 
because the latter obeys an approximate asymptotic relationship for R of the form r''/^ 
(see, for example, Eltayeb & Kumar, 1977). This second transition depends on the value 
of S and will occur at higher values of r and R for lower values of S. There is little 
chance that magneto-inertial convection occurs in the Earth's core, for instance, since S 
is of the order 30000 while the usual estimate for r is 10^^ but it might be relevant for 
understanding of rapidly rotating stars with strong magnetic fields. 

5. Discussion 

A major result of early studies of magnetoconvection in rapidly rotating systems is 
the existence of an absolute minimum critical Rayleigh number (Chandrasekhar, 1961, 
Eltayeb & Kumar, 1977). Imposed magnetic field typically has the effect of increasing the 
azimuthal scale in comparison with non-magnetic convection. When the magnetic field 
is weak the small scale of convection rolls leads to a larger value of the critical Rayleigh 
number as the magnetic field strength decreases. When the magnetic field is strong, it 
inhibits convection and the critical Rayleigh number increases with the increase of the 
strength of the imposed field. Consequently, an optimal state of convection exists at in- 
termediate magnetic field strength. In this paper we consider the weak field regime but 
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a decrease of the critical Raylcigh number with increasing field strength is not observed. 
This is due to the fact that, non-magnetic inertial convection already has a rather large 
azimuthal length scale (Busse & Simitev, 2004) which cannot be increased significantly 
further by the imposed weak magnetic field. 

The support of the Edinburgh Mathematical Society is gratefully acknowledged. 
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